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Abstract: Simple supersymmetric grand unified models based on the gauge group 50(10) 
require -in addition to gauge and matter unification- the unification of t-b-r Yukawa cou- 
plings. Owing to sparticle contributions to fermion self-energy diagrams, the Yukawa uni- 
fication however only occurs for very special values of the soft SUSY breaking parameters. 
We perform a search using a Markov Chain Monte Carlo (MCMC) technique to investigate 
model parameters and sparticle mass spectra which occur in Yukawa-unified SUSY models, 
where we also require the relic density of neutralino dark matter to saturate the WMAP- 
measured abundance. For Yukawa unified models with > 0, the spectrum is characterizd 
by three mass scales: first and second generation scalars in the multi-TeV range, third 
generation scalars in the TeV range, and gauginos in the ~ 100 GeV range. Most solutions 
give far too high a relic abundance of neutralino dark matter. The dark matter discrepancy 
can be rectified by i). allowing for neutralino decay to axino plus photon, ii). imposing 
gaugino mass non-universality or Hi), imposing generational non-universality. In addition, 
the MCMC approach finds a compromise solution where scalar masses are not too heavy, 
and where neutralino annihilation occurs via the light Higgs h resonance. By imposing 
weak scale Higgs soft term boundary conditions, we are also able to generate low /i, tua 
solutions with neutralino annihilation via a light A resonance, though these solutions seem 
to be excluded by CDF/DO measurements of the Bg — > branching fraction. Based 

on the dual requirements of Yukawa coupling unification and dark matter relic density, we 
predict new physics signals at the LHC from pair production of 350-450 GeV gluinos. The 
events are characterized by very high 6-jet multiplicity and a dilepton mass edge around 
m^o — m^o ~ 50-75 GeV. 
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1. Introduction 



Grand unified theories (GUTs) based upon the gauge group SO{W) certainly have to be 
considered as among the most beautiful ideas in particle physics |]l|]. In addition to gauge 
group unification, one also has matter unification of each generation within the 50(10) 
16-dimensional spinorial representation ■0(16). Furthermore, the simplest SO(IO) GUT 
theories also allow for Yukawa coupling unification. The ad-hoc but fortuitous triangle 
anomaly cancellation found in the SM or even SU{5) GUTs is a simple mathematical 
fact in 50(10). The beauty of 50(10) is only enhanced via a marriage to softly broken 
= 1 supersymmetry (SUSY), which stabilizes the gauge hierarchy, and is experimentally 
supported by gauge coupling unification found within the MSSM (provided superpartners 
exist at or around the weak scale). SUSY 50(10) also elegantly addresses the neutrino 
mass problem, since one only has matter unification within the superfield '0(16) provided 
one introduces an additional (SM gauge singlet) superfield N'^ which contains a right- 
handed neutrino state. Under the breaking of 50(10), a superpotential term leading to a 
Majorana mass is allowed for the field; this naturally yields a description of neutrino 
masses in terms of the elegant see-saw mechanism 

Standard GUTs and also SUSYGUTs formulated in 4-d spacetime have fallen into dis- 
repute due to a variety of problems associated with GUT gauge symmetry breaking via the 
Higgs mechanism. These include the doublet-triplet splitting problem, lack of observation 
of proton decay, and the frequently awkward implementation of GUT symmetry breaking 
via at least one large and unwieldy Higgs representation. With the onset of model build- 
ing utilizing extra dimensions, it has been shown to be possible to formulate SUSYGUTs 
in 5 or more spacetime dimensions. Then, the GUT gauge symmetry can be broken via 
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compactification of the extra dimensions on a suitable sub-space, such as an orbifold. In 
these 5-d and 6-d SUSYGUT models, the large GUT scale Higgs representations can be 
dispensed with, the doublet-triplet splitting problem can be solved, and the proton can be 
made longer-lived than current limits or even absolutely stable The extra-dimensional 
SUSYGUT models act as a sort of "proof of principle" of what might be possible in more 
complicated set-ups where the SUSYGUT model might arise from compactification of su- 
perstring models. 

The imminent turn-on of the CERN LHC naturally leads one to ask: How might 
5*0(10) SUSYGUT theories manifest themselves in the environment of an LHC detector? 
Our goal in this paper is to address this question. To do so, our path will be guided by 
the twin requirements of i). Yukawa coupling unification and ii). explaining the measured 
dark matter (DM) abundance of the universe. Our answer we find is that, if DM-allowed, 
Yukawa-unified 50(10) SUSYGUTs are correct, then we expect new physics at LHC to 
consist of gluino pair production events with rrig ~ 350-450 GeV, followed by 3-body gluino 
cascade decays to 6-jet rich final states. In addition, an opposite sign/same-flavor (OS/SF) 
isolated dilepton invariant mass spectrum should have a mass edge at around 50-75 GeV. 
SUSY scalar fields other than the light Higgs h are very heavy, and largely decouple from 
LHC physics. The remainder of this paper details our methodology as to how we come to 
these conclusions. 

In this paper, we assume nature is described by an 50(10) SUSYGUT theory at energy 
scales Q > Mqut ~ 2 x lO^^ GeV. We further assume that the 50(10) SUSYGUT model 
breaks (either via the Higgs mechanism or via compactification of extra dimensions) to the 
MSSM (or MSSM plus right-handed neutrino states) at Q = Mqut- Thus, below Mqut, 
the MSSM is the correct effective field theory which describes nature. We will further 
assume that the superpotential above Mqut is of the form 

/ 9 MeV^iefAio + • • • (1-1) 

so that the three third generation Yukawa couplings /t, /{, and fr are unified at Mqut- It 
is simple in this context to include as well the effect of a third generation neutrino Yukawa 
coupling fy] this effect has been shown to be small, although it can help improve Yukawa 
coupling unification by a few per cent if the neutrino Majorana mass scale is within a few or- 
ders of magnitude of Mqut- Within this ansatz, the GUT scale soft SUSY breaking (SSB) 
terms are constrained by the 50(10) gauge symmetry so that matter scalar SSB terms have 
a common mass mig, Higgs scalar SSB terms have a common mass mio and there is a com- 
mon trilinear soft breaking parameter ^o- As usual, the bilinear soft term B can be traded 
for tan/3, the ratio of Higgs field vevs, while the magnitude of the superpotential Higgs 
mass n is determined in terms of via the electroweak symmetry breaking minimization 
conditions. Here, electroweak symmetry is broken radiatively (REWSB) due to the large 
top quark mass. In order to accomodate REWSB, it is well-known that in Yukawa-unified 
models, the GUT scale Higgs soft masses must be split such that m?u^ < ^'h^ in order to 
fulfill the EWSB minimization conditions; this effectively gives rn'jj^ a head start over rn'jj^ 
in running towards negative values at or around the weak scale. We parametrize the Higgs 
splitting as rrijj- = mfg ^ 2M^. The Higgs mass splitting might originate via a large 
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near-GUT-scale threshold correction arising from the neutrino Yukawa coupling: see the 
Appendix to Ref. for discussion. Thus, the Yukawa unified SUSY model is determined 
by the parameter space 

mi6, mio, Ml, mi/2, ^o, tan/5, sign{fi) (1.2) 

along with the top quark mass. We will take nit = 171 GeV, in accord with recent 
measurements from CDF and DO 

Much previous work has been done on the topic of t~b-T unification in SO{10) SUSYGUTs. 
Early on it was found that t-b-T Yukawa unification could only occur at very high values 
of tan/?|Q. The importance of weak scale MSSM threshold corrections to fermion masses 
was noted by Hall, Rattazzi and Sarid[^. It was also noted that with Yukawa coupling 
unification, in order to obtain an appropriate REWSB, the GUT scale Higgs masses would 
need to be split such that Tnjj^ < fnjj^- perhaps via D-term contributions 1^] to all scalar 
masses (the DT model), or via splitting of only the Higgs soft terms Q (the HS model). 

In Ref. 0, it was found using the Isajet sparticle mass spectrum generator pl| Isasuga 
that Yukawa coupling unification to 5% could be achieved in the MSSM using D-term 
splitting, but only for /i < 0; for fi > 0, the Yukawa coupling unification was much worse, 
of order 30-50%. These parameter space scans allowed miQ values up to only 1.5 TeV, and 
used a GUT scale Yukawa unification quantity 

max{ft, fb, fr) „^ 

— ■ (f f f \ ^ (.1-"^) 

so that e.g. R = 1.1 would correspond to 10% Yukawa unification. The /i < Yukawa 
unification solutions were examined in more detail in Ref. |T^, where dark matter allowed 
solutions were found, and the neutralino j4-annihilation funnel was displayed for the first 
time. 

With the announcement from BNL experiment E-821 that there was a 3(T deviation 
from SM predictions on the muon anomalous magnetic moment a^ = [g — 2)^/2, attention 
shifted back to ^ > solutions. Ref. |l^, using the DT model with parameter space scans 
of mi6 up to 2 TeV, found Yukawa-unified solutions with R '-^ 1.3 but only for special 
choices of GUT scale boundary conditions: 

^0 ~ -2mi6, mio ~ 1.2mi6, (1-4) 

with mii2 <C TTT-ig and tan/5 ~ 50. In fact, these boundary conditions had been found 
earlier by Bagger et al. |l^ in the context of models with a radiatively driven inverted scalar 
mass hierarchy (RIMH) , wherein RG running of multi-TeV GUT scale scalar masses caused 
third generation masses to be driven to weak scale values, while first/second generation soft 
terms remained in the multi-TeV regime. These models, which required Yukawa coupling 
unification, were designed to maintain low fine-tuning by having light third generation 
scalars, while solving the SUSY flavor and CP problems via multi-TeV first and second 



generation scalars. A realistic implementation of these models in Ref. |15] using 2-loop 



RGEs and requiring REWSB found that an inverted hierarchy could be generated, but 
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only to a lesser extent than that envisioned in Ref. [14|, which didn't implement EWSB or 
calculate an actual physical mass spectrum. 

Simultaneously with Ref. [|^, Blazek, Dermisek and Raby (BDR) published results 
showing Yukawa-unified solutions using the HS model solution |Q] . Their results also found 
valid solutions using the Bagger et al. boundary conditions. BDR used a top-down method 
beginning with actual Yukawa unification at Mqut, and implemented 2-loop gauge and 
Yukawa running but 1-loop soft term running. They extracted physical soft terms at scale 
Q = Mz, and minimized a two-Higgs doublet scalar potential to achieve REWSB, also at 
scale Mz- Each run generated a numerical value for third generation t, h and r masses and 
other electroweak and QCD observables. A fit was performed to select those solutions 
which best matched the meaured weak scale fermion masses and other parameters. BDR 
scanned mig values up to 2 TeV, and found best fit results with niA ~ 100 GeV and 
H ~ 100-200 GeV, in contrast to Ref. [|^], where solutions with valid EWSB could only 
be found if /i ~ mA ~ tnii ~ ^ TeV.^ 

In a long follow-up study using Isajet, Auto et al. |18] found that Yukawa-unified 
solutions good to less than a few percent could be found in the /.i > case using the 
HS model of BDR, but only for very large values of mig ^ 5-10 TeV and low values of 
777-1/2 ^ 100 GeV, again using Bagger et al. boundary conditions. Yukawa unification in 
the DT model was at best good to 10% (for this reason, in the present paper we will focus 
only on the HS model). The spectra were characterized by three mass scales: 

1. ~ 5-15 TeV first and second generation scalars, 

2. ~ 1 TeV third generations scalars, /i term and and 

3. chargino masses 777 -± ~ 100-200 GeV and gluino masses 7775 ~ 350-450 GeV. 



These Yukawa-unified solutions- owing to very large values of scalar masses, mA and 
/7- predicted 
result [ll9i of 



/7- predicted dark matter relic density values ^^oh far beyond the WMAP-measured 



Odm/j' = 0.11ll°:[!Ji (2a). (1.5) 

Meanwhile, the spectra generated using the BDR program could easily generate Cl^oh^ 
values close to 0.1 since their allowed /U and 777^ values were far lower, so that mixed 
higgsino dark matter or j4-funnel annihilation solutions could easily be found. In follow-up 



papers to the BDR program [20, 21 1, the neutralino relic density and branching fraction 
Bg — > were evaluated. To avoid constraints on BF{Bs — > from the CDF 

collaboration, the best fit values of tjiiq and 777^ have been steadily increasing, so that 
the latest papers have tua ~ 500 GeV and ttiiq ~ 3 TeV, while fi can still be of order 



100 GeV [21|. In Ref. |22], attempts were made to reconcile the Isajet Yukawa-unified 
solutions with the dark matter relic density. The two solutions advocated were 1. lowering 
GUT scale first/second generation scalars relative to the third, to gain neutralino-squark 



""^A paper by Tobe and Wells [p^ (TW) appeared after Ref. Q. While TW calculate no actual spectra or 
address EWSB, they do adopt a semi-model independent approach which favors t — 6 — r Yukawa coupling 
unification if scalar masses are in the multi-TeV regime while gauginos are as light as possible. 
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or neutralino-slepton co-annihilation solutions [23|, or 2. increasing the GUT scale gaugino 
mass Ml, so the relic density could be lowered by bino-wino co-annihilation |24|. 

In this paper, we report on a new examination of SUSY mass spectra constrained by 
Yukawa coupling unification. We are motivated by the following rationale. 

• We use the latest measured value of the top quark mass- rrit = 171 GeV |||]- whereas 
earlier analyses used mt = 175 GeV or even higher values. 

• We generate SUSY mass spectra using Isajet 7.75, which includes several upgrades 
to the SUSY spectrum calculation from versions used in earlier analyses [p5|. 



• We use a Markov Chain Monte Carlo (MCMC) technique, which is much more effi- 
cient at searching through multi-dimensional parameter spaces. 

• We adopt an alternative approach to implementing Higgs sector SSB boundary con- 
ditions, which simultaneously imposes GUT-scale Higgs (GSH) and weak-scale Higgs 
(WSH) soft term boundary conditions. This allows us to explore solutions akin to 
those found by BDR with low values of rriA and /x. These solutions are difficult to 
generate using a purely top-down approach, due to a high degree of fine-tuning in 
the EWSB sector. 

In Sections |2| and ^, we present results from our calculations. These include: 



1. A variety of solutions consistent with Ref. p8f| with very large values of mie ~ 5-15 
TeV, and very light gauginos. These solutions always have too large a dark matter 
relic density. We propose they may still be valid by invoking an unstable lightest 
neutralino which decays to an axino lightest SUSY particle (LSP). 

2. New solutions are found with niiQ ~ 3 TeV and mi/2 ~ 100 GeV with Yukawa 
unification to ~ 8%, but with a valid relic density due to neutralino annihilation 
through the light Higgs h resonance. 

3. Using simultaneous GSH and WSH boundary conditions, we are able to generate 
BDR-like solutions using Isajet with mig ~ 3-7 TeV but with mi/2 ~ 100-400 GeV 
with good relic density due to neutralino annihilation through the pseudoscalar Higgs 
A resonance. The ^-resonance annihilation solutions have excellent Yukawa unifica- 
tion, but with niA ~ 150-230 GeV and tan (3 ~ 50, are all excluded by measurements 
of the Bg — > fJ.'^/J-'' branching fraction. Solutions with higher tjia and low /i with 
mixed higgsino dark matter also can be found, but tend to have Yukawa unification 
R > 1.2. 

4. We also re-examine the first class of solutions, but solving the DM problem via 
generational non-universality or gaugino mass non- universality, as in Ref. [ p^ . 

In Sec. 1^, we discuss a Table of benchmark scenarios for each of these cases. The cases 
are suitable for collider event generation. We also discuss implications of each of these 
scenarios for collider searches at the CERN LHC. 
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Our preferred solutions- numbers 1. and 2. above- lead to spectra with first/second 
generation scalar masses in the 3 — 15 TeV range, while gluino masses are in the 350 — 450 
GeV range. (Solution number 3 is likely ruled out by the measured limit on BF{Bs — > 
fi^fj,^) and solution 4 gives up either gaugino mass unification or generational unification.) 
Thus, we predict that DM-allowed Yukawa-unified 50(10) SUSY GUT models will lead to 
collider events at the CERN LHC typified by nearly pure gluino pair production followed 
by cascade decays g tbxi, 66x2 ^^'^ ^bxi- The gluino pair events will occur with 
cross sections in the ~ 10^ fb range, and will be very rich in 6-jets. There should exist a 
characteristic OS/SF dilepton mass edge at m^o — m^o ~ 50-75 GeV. 

A summary and conclusions are presented in Sec. ^. 



2. Random scan in HS model 

In this section, we explore the parameter space Eq. (|1.2D for Yukawa- unified solutions by 



means of a random scan. We wish to first check and update results presented in Ref. |18|, 
using the latest Isajet version and a top quark mass of nit = 171 GeV in accord with 
recent measurements from the Fermilab Tevatron Q. The degree of Yukawa unification, 
R, is defined in Eq. ( |1.3D , so that e.g. a value of R = 1.1 corresponds to 10% Yukawa 
unification. 

For our calculations, we adopt the Isajet 7.75 ||ll|, ^] SUSY spectrum generator Isas- 
ugra. Isasugra begins the calculation of the sparticle mass spectrum with input DR gauge 
couplings and fb, fr Yukawa couplings at the scale Q = Mz {ft running begins at Q = rrit) 
and evolves the 6 couplings up in energy to scale Q = Mqut (defined as the value Q where 
Qi = 92) using two-loop RGEs.^ At Q = Mgut, the SSB boundary conditions are input. 



and the set of 26 coupled two-loop MSSM RGEs are used to evolve couplings and SSB 
terms back down in scale to Q = Mz- Full two- loop MSSM RGEs are used for soft term 
evolution, while the gauge and Yukawa coupling evolution includes threshold effects in the 
one- loop beta- functions, so the gauge and Yukawa couplings transition smooothly from the 
MSSM to SM effective theories as different mass thresholds are passed. In Isajet 7.75, the 
SSB terms of sparticles which mix are frozen out at the scale Q = MsusY = y/'^^WJW^^ 
while non-mixing SSB terms are frozen out at their own mass scale |^^. The scalar po- 
tential is minimized using the RG-improved one-loop MSSM effective potential evaluated 
at an optimized scale Q = MsusY which accounts for leading two- loop effects Once 
the tree- level sparticle mass spectrum is computed, full one-loop radiative corrections are 
caculated for all sparticle and Higgs boson masses, including complete one-loop weak scale 
threshold corrections for the top, bottom and tau masses at scale Q = MsuSY- These 
fermion self-energy terms are critical to evaluating whether or not Yukawa couplings do 
indeed unify. Since the GUT scale Yukawa couplings are modified by the threshold cor- 
rections, the Isajet RGE solution must be imposed iteratively with successive up-down 
running until a convergent solution for the spectrum is found. For most of parameter 



^As inputs, we take the top quark pole mass mt = 171 GeV. We also take m^^{Mz) = 2.83 GeV]26[| 
and m?^{Mz) = 1.7463 GeV. The paper Ref. jisj addresses consequences of varying the values of rat and 
rub. 
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space, there is very good agreement between Isajet and the other pubhc spectrum codes 
SoftSusy, SuSpect and SPheno, although at the edges of parameter space agreement be- 
tween the four codes typically diminishes p9t| . 

We first adopt a wide parameter range scan, and then once the best Yukawa-unified 
regions are found, we adopt a narrow scan to try to hone in on the best unified solutions. 
The parameter range we adopt for the wide (narrow) scan is 



(2.1) 



mi6 : 0-20TeV (l-20TeV), 

mio/mie: 0-1.5 (0.8-1.4), 

mi/2 : 0-5 TeV (0-lTeV), 

Ao/miQ -3-3 (-2.5-1.9), 

MD/mie: 0-0.8 (0.25-0.8), 

tan/3: 40-60 (46- 53). 

For the random scan, we evaluate Q^oh'^, BF{h — > 57), Aa^ and BF{Bs — > /U"*"^") using 
Isatools (a sub-package of Isajet). We plot only solutions for which m-± > 103.5 GeV, in 
accord with LEP2 searches, and for the moment implement no other constraints, such as 
relic density, Higgs mass, etc.. 

2.1 Random scan results 

Our first results are shown in Fig. |l], where we show points from the wide scan (dark blue) 
and points from the narrow scan (light blue) in the parameter versus R plane. From frame 
a), we see that Yukawa unification to better than 30% {R < 1.3) cannot be achieved for 
mi6 < 1 TeV, while Yukawa coupling unification becomes much more likely at multi-TeV 
values of m-ig. Frame b) shows that Yukawa-unified models prefer niiQ I — 1.3mi6, while 
frame c) shows that a positive value of ~ (0.25 — 0.5) mig- which yields rnj^^ < rnjj^- 
is preferred. In frame d), we see that the best Yukawa-unified solutions are found for the 
lowest possible values of mi/2- We note here that- using 1-loop RGEs along with the 
LEP2 constraint m-± > 103.5 GeV- one would expect from models with gaugino mass 
unification that since ^tt-^* ~ M2{weak) ~ 0.8mi/2 that we would have mx/2 ^ 125 GeV 
always. However, the very large values of mig we probe alter the simple 1-loop gaugino 
mass unification condition (that — = — = — ) via 2-loop RGE effects. Thus, values of 
mi/2 much lower than ~ 125 GeV are possible if mig is large. 

In frame e), we see a sharp dependence that Yukawa- unified solutions can only be 
obtained for ^0 ~ — 2mi6, while frame f) shows that tan/3 must indeed be large: in the 
range ~ 47 — 53. Bagger et al. had shown in Ref. [^] that a radiatively-driven inverted 
scalar mass hierarchy with m{third generation) <C m{fir st / second generation) could be 
derived provided one starts with unified Yukawa couplings, the boundary conditions 

4ml^ = 2mjo = Al (2.2) 

and one neglects the effect of gaugino masses. Our results in Fig. [I| show the inverse effect: 
that Yukawa coupling unification can only be achieved if one imposes the boundary condi- 
tions ( |2.2[ ) along with mig ^ fni/2- This result holds only in our numerical calculations for 
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/X > and Aq < and of course m'^j^ < m^u^ ■ The results shown in Fig. |T| also verify that 



the results obtained in Ref. |18| still hold, even with updated spectra code and a lower 
value of mt = 171 GeV. 

In Fig. ^, we show various -ino masses^ versus R as generated from our random scan. 
In frame a), we see that- owing to the preference of Yukawa-unified solutions to have 
as small as possible, the chargino mass m-± is preferred to be quite light, as close to the 
LEP2 limit as possible, with m-± ~ 100-200 GeV. Likewise, in frame b), the gluino mass 
should be relatively light, with rrig ~ 350-500 GeV. The lightest neutralino Xi mass is 
shown in frame c), and is preferred in the range m^o ~ 50-100 GeV. Meanwhile, the mass 
difference mj^o — m^o is shown in frame d), and is also in the range ~ 50-100 GeV. This 
latter quantity is important because if m^o — m^o < Mz, two body spoiler decay modes 
such as X2 ~^ Xi^ will be kinematically closed, and the three body decays X2 ~^ Xi^^ 
{i = e or 11) should occur at a sufficiently large rate at the LHC that an edge should be 
visible in the m{U) invariant mass distribution at m^o — m^o |^]. This measureable mass 
edge can serve as the starting point for sparticle mass reconstruction in SUSY particle 
cascade decay events at the LHC [Q. Thus, in Yukawa-unified models, this mass edge is 
highly likely to be visible. 

In Fig. ^, we show the expected masses of a) UL-squark, h) the fi-squark, c), the pseu- 
doscalar Higgs boson A and d) the superpotential Higgs parameter //. Frame a) shows that 
Yukawa-unified solutions prefer first/second generation squarks and sleptons with masses 
in the 5-20 TeV range- far higher than values typically examined in phenomenological 
SUSY studies! The top squark mass and the A, H and Higgs bosons tends to be 
somewhat lighter: in the 2-8 TeV range. Finally, frame d) shows that the ^ parameter- 
which is derived from the EWSB minimization conditions- tends also to be in the 5-15 
TeV range. Thus, using a top-down approach to search for Yukawa-unified solutions in the 
HS model, we find that /i ^ Mi, M2, so that the lighter charginos and neutralinos should 
be gaugino-like, and quite light, while the heavier charginos and neutralino will be in the 
muti-TeV range, and nearly pure higgsino-like states. In particular, the lightest SUSY 
particle (LSP)- the neutralino is nearly pure bino-like. 

In Fig. ^, we plot R vs. ^^oh'^ for LEP2 allowed points from our random scan. It is 
clear that ~ 1 points predict an extremely large value of O^o/i^ of 30-30,000. On the 

Al 

other hand, if we require consistency with the WMAP-measured value of ^j^oh'^ ~ 0.1, then 
we generate Yukawa-unified solutions to 40% unification with the random scan. This plot 
underscores the difficulty of finding sparticle mass spectra solutions which are compatible 
with both the measured dark matter abundance and t-b~T Yukawa coupling unification. 



^We collectively refer to the set of all gluinos, charginos and neutralinos as -inos. 
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Figure 1: Plot of R versus various input parameters for a wide (dark blue) and narrow (light blue) 



random scan over the parameter ranges listed in Eq. (2.1). We take fi > and rut = 171 GeV. 
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Figure 2: Plot of R versus various sparticle masses for a random scan over tire parameter range 
listed in Eq. (O). We take and mt = 171 GeV. 
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Figure 3: Plot of R versus various sparticle masses for a random scan over the parameter range 
listed in Eq. (|^). We take ^ > and rrit = 171 GeV. 
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2.2 Three proposals to reconcile Yukawa- unified models with dark matter relic 
density 

2.2.1 Dark matter solution via neutralino decay to axino 

We see from Fig. ^ that models generated from the random scan with i? ~ 1.0 aU have 
n^oh'^ ~ 30 - 30,000: far beyond the WMAP-measured result of ^cDuh^ ~ 0.1. One 
possible solution to reconcile the predicted and measured dark matter density is to assume 
that the lightest neutralino is in fact not the LSP, but is unstable. Some alternative 
LSP candidates consist of the gravitino G or the axino a. In gravity-mediated SUSY 
breaking models, the gravitino mass 772-3/2 arises due to the superHiggs mechanism, and is 
expected to set the scale for all the soft SUSY breaking terms. Usually it is assumed the 
gravitino is heavier than the lightest neutralino rn-3/2 > ^j^Oj in which case the gravitino 
essentially decouples from phenomenology. However, if 7113/2 < f^yfi-, then the x? becomes 
unstable and can decay via modes such as Xi ~^ iG. The Xi lifetime is expected to be 
very long- of order 10^ — lO^'^ sec- so the neutralino still escapes detection at collider 
experiments, but is susceptible to constraints from Big Bang nucleosynthesis (BBN) and 
CMB anisotropics ||3^. The relic density of gravitinos is expected to be simply Vl^ = 

^^O^o/i^, since the gravitinos "inherit" the thermally produced neutralino relic number 

^1 

density. Thus, a scenario with a G super WIMP as LSP in SUGRA-type models can reduce 
the relic density by typically factors of a few- which is not enough in the case of Yukawa- 
unified models, where relic density suppression factors of 10^ — 10^ are needed. 

A better option occurs if we hypothesize an axino a LSP. If indeed there is a Peccei- 
Quinn solution to the strong CP problem, then one expects the existence of axions, typi- 
cally with mass below the eV scale. While axions can themselves form cold dark matter, 
it is also easily possible that they contribute little to the CDM relic density. However, in 
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models with SUSY and axions, then the axion is just one element of an axion superfield, 
the superpartner of the axion being a spin-^ axino a. The axino mass can be far different 
from the typical soft SUSY breaking scale, and the range rria ~ eV-GeV is allowed. 

Axinos can be produced in the early universe both thermally or non-thermally from 



NLSP decay. From the latter source, we expect roughly |33| 



n- ^ H^Q oh"^. (2.3) 
Thus, for Q^o ~ 10'^ and with m^o ~ 50 GeV as in Yukawa-unified models, an axino mass 

Xl Xl 

of rria ^ 5 MeV is required 

In this mass range, the axinos from Xi decay are expected to give a hot /warm com- 



ponent to the dark matter |34|. However, thermally produced axinos in this mass range 
could yield the required cold dark matter. Thus, if an unstable neutralino decay Xi ~^ 07 
is to reconcile Yukawa-unified models with the relic density, then we would expect the 
dark matter to be predominantly cold axinos produced thermally, but with a re-heat tem- 
perature Tr < Tf, where Ty is the temperature where axinos decouple from the thermal 
plasma in the early universe. This scenario admits a dark matter abundance that can be in 
accord with WMAP measurements, and would be primarily CDM, but with a warm dark 
matter component arising non-thermally from Xi decays. For a bino-like neutralino, as in 
Yukawa-unified models, the Xi lifetime is given by [^5| 

2 1 f fa/N \V50GevV , 
,.3.3xl0-sec^(^) (^—j , ,2.4) 

where the model-dependent constant CaYY is of order 1, is the Peccei-Quinn breaking 
scale, and is a model dependent factor (A^ = 1(6) for the KSVZ (DFSZ) axion model). 
Thus, for reasonable choices of model parameters, we expect the neutralino lifetime to be 
of order 3 x 10~^ sec. This is short enough so that photon injection into the early universe 
from Xl — * 07 decay occurs before nucleosynthesis, thus avoiding constraints from BBN. 

For illustration, we adopt a point A listed in Table |l] of Yukawa-unified benchmark 
models. The point has mie = 9202.9 GeV, mio = 10966.1 GeV, Md = 3504.4 GeV, 
mi/2 = 62.5 GeV, Aq = -19964.5 GeV, tan/3 = 49.1 GeV with ^ > and = 171 GeV. 
It has m^o = 55.6 GeV and ^^^oh'^ = 423 (IsaReD result). Thus, Xi ~^ '^7 with rria ^ 1 
MeV would allow for a mixed warm/cold axino dark matter solution to the problem of relic 
density in Yukawa-unified models. 

2.2.2 Dark matter solution via non-universal gaugino masses 

An alternative solution to reconciling the dark matter abundance with Yukawa-unified 
models is to consider the possibility of non-universal gaugino masses. If we adopt any of 
the Yukawa unified models from the random scan and vary the SU{2) {SU{3)) gaugino 
masses M2 {M3), then the Yukawa coupling unification will be destroyed via the effect of 
iixf {go) loops. However, if Mi is varied, Yukawa coupling unification is preserved since 
contributions to fermion masses from loops containing are small. 
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Figure 5: Plot of variation in fl^oh"^ versus non-universal GUT scale gaugino mass Mi for bench- 
mark point A in Table |l|. 



By raising the GUT scale value of Mi to values higher than mi/2, the weak scale value 
of Ml is also increased. If Mi is increased enough, then m^o (which is nearly equal to Mi 
since Xi is largely bino-like) becomes close to ^--i. When this happens, the Xi becomes 



more wino-like, with an increased annihilation cross section to WW pairs if m^o > M\y |]36| . 
In our case, usually m^o < M^y. Then raising Mi still lowers the relic density, but now 

Al 



via bino-wino co-annihilation (BWCA) p4[ . 

In Fig. I we show the variation in ^^oh versus Mi{Mgijt) for benchmark point A in 
Table ||. The location of Mi for point A is marked by the arrow. The double dips at low Mi 
are due to neutralino annihilation through the Z and h poles. Once Mi [Mgut) is increased 
to ~ 195 GeV, then we reach a relic density in accord with WMAP measurements. Since 
m-± ~ m^o, and m-± ~ m^o, the Vn — ¥? mass gap is small, of order 10-20 GeV. We list 

Xl Al Xi Al 

the raised Mi{Mgut) = 195 GeV point as point B in benchmark Table |l[ 

2.2.3 Dark matter solution via generational non-universality 

Another possibility for reconciling the neutralino relic density with the measured value is 
to lower the first/second generation scalar masses mi6(l,2), while keeping mi6(3) fixed 
at miQ. The Bagger et al. inverted hierarchy solution depends only on third generation 
scalar masses, while the effects of the first two generations decouple. Ordinarily, solutions 
with mi6(l, 2) = mi6(3) are taken to enforce the super-GIM mechanism for suppression of 
flavor changing neutral current (FCNC) processes. Limits from FCNCs mainly require near 
degeneracy between the first two generations, while limits on third generation universality 
are much less severe |37]. Lowering mi6(l,2) works to lower the relic density because of 
the large S term in the scalar mass RGEs: 

S = n4j^- mjj^ + Tr [m^ - mi - 2m^ + + m|] . (2.5) 
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Figure 6: Plot of variation in fl^oh'^ versus non- universal GUT scale first/second generation scalar 
mass TOi6(l,2) for benchmark point C in Table |^. 



In models with universaUty, Uke mSUGRA, 5 = to one-loop at all energy scales; in models 
with non-universal Higgs scalars, like the HS model, this term can be large and have a major 
influence on scalar mass running. The large S term helps suppress right-squark masses. If 
mi6(l,2) is taken light enough, then m^^^ ~ — "^x?' neutralino-pair annihilation 
into quarks and neutralino-squark co-annihilation can act to reduce the relic density. 

In Fig. ^ we show the variation in ^^oh'^ versus mi6(l,2) where we take mi6(3) = 
5018.8 GeV, mi/2 = 160 GeV, = -10624.2 GeV, tan/5 = 47.8 and /i > 0. When 
mi6(l, 2) is lowered to 603.8 GeV, then mn^ ~ = 98.3 GeV, and we have neutralino an- 
nihilation via light i-channel squark exchange and also neutralino-squark co-annihilation.'* 
IsaReD and Micromegas give Q.^oh'^ ~ 0.1 at this point, which we adopt as benchmark 
point C in Table ||. The two light squarks are just at the limit of LEP2 exclusion. They 
may possibly be excludable by Tevatron analyses, but the squark-neutralino mass gap is 
quite small, so the energy release from uji — > uxi is low. So far, no such study has been 
made, and so the possibility cannot yet be definitively excluded. 



*A bug fix is needed in the Isajet 7.75 IsaReD subroutine in order to obtain the correct relic density; we 
thank A. Pukhov for pointing this out. 
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3. Markov Chain Monte Carlo analysis 



Next we adopt an improved scanning method based on the Markov Chain Monte Carlo 
(MCMC) technique to search more efficiently for parameter space regions of good Yukawa 



unification and WMAP-compatible DM relic density. A Markov Chain [38| is a discrete- 
time, random process having the Markov property, which is defined such that given the 
present state, the future state only depends on the present state, but not on the past states. 
That is: 

P(X*+^ = = x\ = x^) = = = xt). (3.1) 

An MCMC constructs a Markov Chain through sampling from a parameter space with 
the help of a specified algorithm. In this study, we have applied the Metropolis-Hastings 
algorithm ||3^ , which generates a candidate state x'^ from the present state x* using a 
proposal density x'^). The candidate state is accepted to be the next state x*^^ if the 
ratio 

^ Pix')Qix';x') 
^ P(x*)Q(x'=;x*)' ^ ' 

(where P{x) is the probability calculated for the state x) is greater than a uniform random 
number a = U{0, 1). If the candidate is not accepted, the present state is retained and a 
new candidate state is generated. For the proposal density we use a Gaussian distribution 
that is centered at x* and has a width a. This simplifies the p ratio to P[x^) / P{x^). 

Once taking off from a starting point, Markov chains are aimed to converge at a target 
distribution P{x) around a point with the highest probability. The time needed for a 
Markov chain to converge depends on the width of the Gaussian distribution used as the 
proposal density. This width can be adjusted during the run to achieve a more efficient 
convergence. 

In our search in the 50(10) parameter space, we assume flat priors and we approximate 
the likelihood of a state to be e~^^^^\ We define the for R as 

xi- (3.3) 

where Runi fication = 1 and aR is the discrepancy we allow from absolute Yukawa unification, 
which, in this case we take to be 0.05. On the other hand, for il/i^ we define 



1, (0.094 < nh"^ < 0.136) 

{nh"^ < 0.094 or nh"^ > 0.136) 



where Qh'^f,^^ = 0.115 is the mean value of the range 0.094 < nh^ < 0.136 proposed 



in |40|, and cTQh'^ = 0.021. This way, the MCMC primarily searches for regions of Yukawa- 
unifications, and within these regions for solutions with a good relic density. 

For each search, we select a set of ~ 10 starting points in order to ensure a more thor- 
ough investigation of the parameter space. Then we run the MCMC, aiming to maximize 
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the likelihood of either R alone, or R and O/i^ simultaneously. For the case of simultaneous 
maximization, we compute the p ratios for R and U.h? individually, requiring both p/j > a 
and Pn/i2 > a separately. We do not strictly seek convergence to an absolute maximal 
likelihood, but we rather use the MCMC as a tool to reach compatible regions and to 
investigate the amount of their extension in the SO (10) parameter space. 



3.1 HS model: neutralino annihilation via h resonance 

We begin our MCMC scans by selecting 10 starting points "pseudorandomly" -that is, 
selecting them from different mig regions to cover a wider range of the parameter space- 
and imposing some loose limits (defined by our previous works and random scans) on the 
rest of their parameters to achieve a more efficient convergence. Our initial scan is directed 
to look for points only with R as close to 1.0 as possible by maximizing solely the likelihood 
of R. Based on the results of the first MCMC scan, we then pick a new set of 10 starting 
points with low R and also low il^o/i^, and direct the second scan to look for points with 

Al 

both R = 1.0 and < 0.136 by maximizing the likelihoods of R and fi/i^ simultaneously. 

Al 

For MCMC scans, the code is interfaced to the micrOMEGAs [^] package to evaluate the 
relic density and low-energy constraints. 

Figure shows the Yukawa-unified region found by the MCMC results as a projection 
in the plane of mig versus m-ig. The light-blue dots are points which have R < 1.1, while 
dark blue dots have R < 1.05. In addition, we show in orange (red) the points which satisfy 
R < 1.1 (1.05) and ^^oh^ < 0.136. The points with low R are narrowly correlated along 
the line mio — 1.2miQ. While the low R points range over mig values from 3 to over 12 
TeV (in agreement with the results from the random scans) the MCMC has also identified 
a range of points with both ~ 1 and fl^oh'^ < 0.136, but only for mie values of about 3-4 

Al 

TeV! 

Fig. Q shows the MCMC scan results in the mie vs. Ao/miQ plane. Again, we see that 
points with low R populate the region with Aq ~ (2-2.1)mi6 over a wide range of mig 
values. The plot includes the Q^oh'^ < 0.136 points around miQ ~ 3-4 TeV. 

Al 

In Fig. ^, we show MCMC results in the miQ vs. mi/2 plane. Here, we see the very 
lowest R points select out the lowest possible ini/2 values allowed for a given value of 
mi6, and that the minimum value allowed steadily decreases with increasing mig- 

the boundary being determined by the LEP2 limit on chargino masses. The points with a 
"good" relic density are clustered around mi/2 ~ 100 GeV. 

We also show in Fig. |lO| the individual GUT-scale values of Higgs soft terms tuh^ 
(lower branch) and mn^ (upper branch). This plot displays the required Higgs splitting 
and confirms that mjy^ > mH„- 



In Fig. |11|, we show points with low R in the ruh — 2m^o vs. ruA — 2?7i^o plane. In 



these solutions, tua is usually far greater than 2m^o, indicating the neutralino annihilation 
through the ^-resonance is not the cause of the reduced relic density orange and red 
points. However, the low ^j^o/i^ points all do lie along the nih — 2m^o line, indicating that 
/i-resonance annihilation is the mechanism at work to reduce the relic density in the early 



universe. In Fig. W^L we show R vs. Vl^oh? for the MCMC scan. In this frame, we see that 

I 1 Al 
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the points with high rehc density extend down to -R = 1, while the low relic density points 
reach below R = 1.05, but can reach no lower than R = 1.03. 

We are now in a position to understand this new class of Yukawa-unified, DM-allowed 
solutions: the search for low R pushes mig to very high, multi-TeV values. Meanwhile, in 
order for /i-resonance annihilation to reduce the relic density to the WMAP-allowed range, 
mi6 can't be too large. The region around mie ~ 3-4 TeV offers a compromise between 
these two tendencies: for mie not too large, the dip in relic density due to the /i-resonance 
annihilation is sufficient to bring the relic density into the desired range. But since mig 
can't be too large, the Yukawa unification is limited to a couple of percent at best. This 
new class of solutions was difficult to reach using a random scan, since the /i-resonance is so 
narrow. The necessary value of m^o has to be just right- with 2m^o slightly below mh~ so 
that the thermal averaging of neutralino energies convolutes with the resonant cross section 
with enough strength to give substantial neutralino annihilation in the early universe. 

The 50(10) model parameters leading to low R and good relic density occur only over 
a very narrow range of w,]^/2 ~ 100 GeV and mig ~ 3 TeV. This means the Yukawa- unified 
/i-resonance annihilation points have very specific mass spectra predictions. We show in 
Fig. ^ the rrig vs. m^^ plane for MCMC points in the HS model. Here, we see that the points 
with n^oh"^ < 0.136 ah have rrig ~ 350-450 GeV, while m^^ ~ 350-750 GeV. The large ^ 
parameter combines with gaugino mass unification to predict that 'm^± ~ m^o ~ 100-150 
GeV, while m^o ~ 50-75 GeV. 

Given this very tightly correlated mass spectrum, gluino cascade decay events at the 
LHC will lead to X2 production, followed mainly by X2 ~^ Xibt> decay and also by X2 ~^ Xi^^ 
decay. The OS/SF isolated dilepton mass spectrum will be bounded kinematically by 
m^o — m^o . The mass difference provides an edge in the dilepton mass spectrum which is 
characteristic of these decays, and which is easily measureable. Furthermore, it should be 
correlated with ruh , since all these masses are related by resonance annihilation and theory. 
The value of rrih should be directly measureable at LHC after several years of data taking 
via the bump in the /i — > 77 mass spectrum. Thus, we plot in Fig. 14 our Yukawa-unified 
MCMC points in the m^o — m^o vs. m/j plane. For this scenario to be bourne out, we 

A,2 Xl 

would predict m^o — m^o ~ 52-65 GeV, with a roughly linear correlation with m/^. 

A2 Al 

We adopt point D in Table |l| as being representative of the light Higgs /i-resonance an- 
nihilation compromise solutions. The relic density computed with micrOMEGAs (il^o/i^ = 
0.06) is below the preferred range, while IsaReD gives ^^oh? = 0.1. Yukawa couplings are 
unified at the 9% level. We note here that we could have adopted a solution with even 
better Yukawa coupling unification at the 4-5% level. These solutions tend to give light 
Higgs mass nih ^110 GeV (as can be seen by the red dots in Fig. 14) which are more 
likely to be excluded by LEP2 Higgs search results. 



-17- 




2000 



2000 4000 6000 8000 10000 
mi6 

Figure 7: Plot of MCMC results in the mig vs. toiq plane; the light-blue (dark-blue) points have 
R < 1.1 (1.05), while for the orange (red) points R < 1.1 (1.05) and n^oh'^ < 0.136. 
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Figure 8: Plot of MCMC results in the mig vs. Ao/miQ plane; the light-blue (dark-blue) points 
have R < 1.1 (1.05), while for the orange (red) points R < 1.1 (1.05) and fi^o/i^ < 0.136. 
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Figure 9: Plot of MCMC results in the mig vs. mi/2 plane; the light- blue (dark-blue) points have 
R < 1.1 (1.05), while for the orange (red) points R < 1.1 (1.05) and n^oh"^ < 0.136. 
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Figure 10: Plot of MCMC results in the mie vs. „ plane; the light-blue (dark-blue) points 
have R < 1.1 (1.05), while for the orange (red) points R < 1.1 (1.05) and n^oh'^ < 0.136. 
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Figure 11: Plot of MCMC results in the m^ — 2m^o vs. mA — 2m^o plane; the light-blue (dark-blue) 
points have R < 1.1 (1.05), while for the orange (red) points R < 1.1 (1.05) and ^l^oK^ < 0.136. 
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Figure 12: Plot of MCMC results in the R vs. il-^ah^ plane; the light-blue (dark-blue) points have 
R < 1.1 (1.05), while for the orange (red) points R < 1.1 (1.05) and n^oh'^ < 0.136. 
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Figure 13: Plot of MCMC results in the nig vs. to^^ plane; the light-blue (dark-blue) points have 
R < 1.1 (1.05), while for the orange (red) points R < 1.1 (1.05) and n^oh'^ < 0.136. 
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Figure 14: Plot of MCMC results in the m^o — m^o vs. irih plane; the light-blue (dark-blue) 
points have R < 1.1 (1.05), while for the orange (red) points R < 1.1 (1.05) and il^oh^ < 0.136. 
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3.2 Solutions using weak scale Higgs boundary conditions 

In the analysis put forth by BDR Q, Yukawa- unified solutions are found with low values 
of both ^ and in the 100-200 GeV range, while mig and mio are typically at the 2-3 
TeV scale. We have seen from our results so far that ^ and tua are typically in the TeV 
regime. Some low solutions were generated using Isajet in Table 2 of Ref. but these 
had R ~ 1.25. 

We find here that we can generate small /i and small solutions using Isajet by 
using the pre-programmed non-universal Higgs model (NUHM)^. The approach is to start 
with a set of GSH soft term boundary conditions, and evolve the soft SUSY breaking Higgs 
masses rn?H^ and m|^^ down to the weak scale MsjjSY- At Q = MsjjSY, re-calculate what 
and ttih^ should have been in order to get the input values of rriA and /i, using the 
two electroweak symmetry breaking minimization conditions (in practise, we use 1-loop 
relations) : 

(mir +m'jr +2^2) sin 2/? 
B = ^ — and (3.5) 



2/i 

^2 ^ IIH^ 7^"^''"^ ^ - (3.6) 



(tan2/?- 1) 2" 



Then run back up to the GUT scale using these new WSH boundary conditions. At 
each iteration, the weak scale values of rnjj^ and w-f^^ have to be re-computed so as to 
maintain the input value of fi and m^; in this case, the GUT scale values of rn'jj^ and 'm-'j^^ 
are outputs, instead of inputs. For this class of solutions, both GSH and WSH boundary 
conditions must be used in Isajet. The GSH boundary conditions are needed just to get 
an acceptable EWSB on the first iteration so that a spectrum can be computed, and then 
modified to yield the input values of and fj,. Using default universal GSH soft terms 
will usually fail to give appropriate EWSB on any iteration where Yukawa couplings are 
unified. 

We implement an MCMC scan over the modified parameter space 

"116, "11/2, ^0, tan/3, TUA, /U (3.7) 

(effectively trading the GUT scale inputs m'jj^ and n^'jj^ (or alternatively mio and M-^) for 
weak scale inputs and ^). We begin with 10 starting points selected pseudorandomly 
from different regions of the above parameter space, and implement two MCMC scans on 
them, one searching for points with lowest R values by maximizing the likelihood of R and 
the other for solutions with R = 1 and ^■^oh'^ < 0.136 by maximizing likelihoods of R and 
O/i^ simultaneously. 

Our first results are shown in Fig. |l5|for the mie vs. Ao/mie plane, where we plot points 
with R < 1.1 (1.05) using dark blue (light blue) dots, and solutions with ^j^oh"^ < 0.136 for 
R < 1.1 (1.05) using orange (red) dots. While we again get good Yukawa-unified solutions 
over a wide range of multi-TeV values of mig, this time we pick up additional dark matter 
allowed solutions for tuiq : 3-6 TeV. The solutions again respect the Bagger et al. boundary 
condition Aq ~ —2miQ. 

^This is model line 8 of the Isajet non- universal supergravity models (NUSUG) 
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In Fig. 16, we show the WSH solutions in the mig vs. mi/2 plane. The minimum in 
allowed mi/2 values again decreases with increasing mig. We see that for the WSH class 
of solutions, much larger values of mii2 ranging up to 300 — 500 GeV are DM-allowed. 



In Fig. 17, we plot the WSH solutions in the input parameter niA vs. /i plane. In this 
case, we see that the bulk of the DM-allowed solutions occur at relatively low values of 
TTLA ~ 130-250 GeV. These low rriA solutions were extremely difficult to generate with the 
top-down approach, and indicate that they have a high degree of fine-tuning.^ A scattering 
of DM-allowed dots occur with high mA values. These turn out to be the /i-resonance 



solutions as generated with the GSH boundary conditions in Sec. 3.1 

This is seen more clearly by plotting in the m/i — 2m^o vs. rriA — 2m^o plane in Fig. W& 

K\ K\ I— 

Here we see a narrow strip at rrih — 2m^o = corresponding to /i-resonance annihilation 
solutions, while we also have a wider band of solutions at niA — 2m^o = 0, which indicate 

Al 

neutralino annihilation through the A-resonance. The width of the latter band is due to 
the fact that the A width can be quite wide- typically a few GeV, while the /i-width is 
much narrower, of order 50 MeV. 

The A-resonance solutions occur at tan /3 ~ 50 and relatively low ruA values. This can 
signal dangerously high branching fractions for A*^/^^ decay [^] since the branching 

fraction goes like tan^ . We plot the BF{Bs fi'^fi^) vs. m/i in Fig. |19[ The 
recent experimental limit from the CDF collaboration is that BF{Bs — > fi~^^~) < 5.8 x 
10-8 11]. Thus, the entire band of A -resonance annihilation solutions becomes excluded! 
The smattering of DM-allowed dots below the CDF limit all occur with DM annihilation 
via the /i-resonance. 

In case these A-resonance solutions are somehow allowed- say by additional flavor- 
violating soft terms- we plot the solutions in the rrig vs. m^o plane in Fig. Here, we 
see a much larger range of rrig and m^o values are DM-allowed than in the GSH solutions, 
with nig extending up to 1500 GeV. 

If the DM-allowed GSH solutions are able to avoid the BF{Bs fx^fj,^) constraint, 
then the values of niA and m^o — m^o will be correlated, and the latter quantity will be 
measureable if the mass gap m^o — m^o < Mz- The predicted correlation is shown in 
Fig. In this case, the mass gap m^o — m^o runs far beyond Mz- 

We present a benchmark point with low niA and ~ 1 as point E in Table ||. While 
the point is DM-allowed, it also violates the CDF bound on BF{Bs — > 

At this point, it is useful to compare the Isajet SUSY spectral solutions to those 
generated by Dermisek et al. in Ref. pO] and [21|. In Fig. l2^, we plot the Isajet 7.75 



solutions in the mi/2 vs. /x plane for tjiiq = 3 TeV, mio/mig = 1.3, Ao/niiQ = —1.85, 
tan /? = 50.9 and Amjj- = 0.14, with uia = 500 GeV: i.e. corresponding closely to Fig. 



1 of IQ. We plot contours of R from 1.15 to 1.3. Also, the green-shaded regions give 



the WMAP-measured relic density, while white-shaded regions give n^oh^ < 0.095, and 

Al 

pink-shaded regions give ^^oh'^ > 0.13 (as in Dermisek et al). The LEP2 constraint on 
m-± is indicated by the solid contour at low mi/2 and low fj,. We see qualitatively the same 

Al ' 



®The TW paper (Ref. jl^]) remarks that there must be considerable fine-tuning as well to reconcile 
BF{b — > S7) with Yukawa unification and the dark matter relic abundance. 



-23- 



shape to the DM-allowed regions as generated by Dcrmisek et al: the thick green regions 
are DM-allowed either by 74-resonance annihilation at large n, or by mixed higgsino DM 
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Figure 15: Plot of MCMC results using WSH boundary conditions in the mig vs. Ao/mie plane; 
the light-blue (dark-blue) points have R < 1.1 (1.05), while for the orange (red) points i? < 1.1 (1.05) 
and n^oh'^ < 0.136. 
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Figure 16: Plot of MCMC results using WSH boundary conditions in the mig vs. mi/2 plane; the 
light-blue (dark-blue) points have B < 1.1 (1.05), while for the orange (red) points R < 1.1 (1.05) 
and n^oh'^ < 0.136. 
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Figure 17: Plot of MCMC results using WSH boundary conditions in the uia vs. /i plane; the 
light-blue (dark-blue) points have R < 1.1 (1.05), while for the orange (red) points R < 1.1 (1.05) 
and n^oh'^ < 0.136. 




Figure 18: Plot of MCMC results using WSH boundary conditions in the mh — 2m^o vs. ■mA — '2m^o 
plane; the light-blue (dark-blue) points have R < 1.1 (1.05), while for the orange (red) points 
R < 1.1 (1.05) and fl^oh^ < 0.136. 
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Figure 19: Plot of MCMC results using WSH boundary conditions in the vs. BF{Bs /^^/i ) 
plane; the light-blue (dark-blue) points have R < 1.1 (1.05), while for the orange (red) points 
R < 1.1 (1.05) and ^o/i^ < 0.136. 
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Figure 20: Plot of MCMC results using WSH boundary conditions in the nig vs. m^o plane; the 
light-blue (dark-blue) points have R < 1.1 (1.05), while for the orange (red) points R < 1.1 (1.05) 



and n^o/i2 < 0.136. 
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Figure 21: Plot of MCMC results using WSH boundary conditions in the m^o — m^o vs. uia 
plane; the light-blue (dark-blue) points have R < 1.1 (1.05), while for the orange (red) points 
R < 1.1 (1.05) and fl^oh'^ < 0.136. 




Figure 22: Contours of R and DM-allowed regions in the mi/2 vs. /z parameter space for mie = 3 
TeV, mio/mi6 = 1.3, Ao/mie = -1.85, Amn = 0.14, tan/3 = 50.9, ttia = 500 GeV and mt = 173.9 
GeV, as in Dermisek et al., but using Isajet 7.75 for mass spectra generation. 
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annihilation at low /x. There is also a light Higgs /i-resonance solution at ~ 120 GeV. 



A notable feature of Fig. 22 is that over much of the DM-allowed region, the Yukawa 
unification has R > 1.27 As we move to larger n values and lower mi/2 values, the Yukawa 
unification gets better and better. Most of the region with i? < 1.15 is DM- forbidden, save 
for the upper part of the light /i-resonance solution. In fact, now we can see why our 
compromise solution (point D) works and why it is so hard to find using a top-down 
approach: only the very narrow upper tip is both DM-allowed, and has a low R value. 

4. Yukawa-unified benchmark scenarios and LHC signatures 

We have assembled in Table |l| five Yukawa-unified benchmark scenarios that yield the 
correct relic abundance of dark matter in five different ways. With the LHC turn-on being 
imminent, it is fruitful to examine what each of these five scenarios implies for new physics 
signatures. 

At the bottom of Table |l] we list Q^oh'^, BF{b 57), BF{Bs — > jjL'^ IS.a^ and spin- 
independent neutralino-proton direct DM detection cross section ct(xiP)- For the first four 
of these numbers, we list output from IsaReD/Isatools (upper) and micrOMEGAs (lower). 
While the results for the low-energy constraints agree fairly well, there is almost a factor 
of 2 difference in the relic density when the neutralino dominantly annihilates through h 
or A exchange (points A, D, E). This is due to differences in the treatment of the Higgs 
resonance. For example, IsaReD in Isajet 7.75 uses Yukawa couplings evaluated at scale 
Q = yjmii^ annihilation through the A resonance and for evaluation of the heavy 

Higgs widths, while micrOMEGAs uses an effective Lagrangian approach and Q = 2m^o.^ 

Point A 

Point A of Table Q is a generic Yukawa-unified model with first and second generation 
scalar masses ~ 9 TeV, so they essentially decouple from LHC physics. Third generation 
and heavy Higgs scalars have masses at the 2-3 TeV level, while the lightest charginos, 
neutralinos and gluinos all have masses in the range 100-400 GeV. Since Jl^o/i^ ~ 400, we 
postulate that the neutralino Xi is in feet an NLSP, decaying to 07 with a lifetime of order 
0.03 seconds. In this case, the mean decay distance of a Xi '^iii be of order 10^ km. Thus, 
the Xi '^iil still escape the LHC detectors, leading to missing energy signatures (although 
it is conceivable some may decay occassionally within the detector). 



The LHC SUSY events will consist of a hard and soft component |44|. The hard 

component comes from pair production of ~ 400 GeV gluinos. The gluinos decay via 3- 

body modes dominantly via g — > tbxf, bbxi ^'^d especially 66x2 ^^l- The gg production 

^Note that although the general features in Fig. here and Fig. 1 of ji^] are similar, the latter results 
were obtained in a top-down fit to low energy obervables assuming exact Yukawa unification, which is a 
different approach then the one followed here. Moreover, there are several important differences in the 
level of sophistication of the spectrum computations between Ref. |2^, |2^] and the study presented 
here. For instance, Ref. H |l) has only 1-loop RGE running of the SUSY-breaking parameters, takes 
sparticle masses to be running masses at scale Q = Mz; ISAJET 7.75 applies full 2-loop running plus 1-loop 
threshold corrections. 

*A complete discussion of the details of the calculations in the two programs is beyond the scope of this 
paper; we refer the interested reader to the respective manuals. 
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cross section is of order 10^ fb at LHC, so we might expect 10^ gluino pair events per 100 
fb~^ of integrated luminosity. After cascade decays, we expect an assortment of events 
with high jet and 6-jet multiphcity, plus an assortment of isolated leptons. The X2 ~^ Xi^e 
branching fraction is at 2.2% , which should be enough to reconstruct the dilepton mass 
edge at m^o — m^o ~ 73 GeV. Correct pairing of 6-jets and/or 6-jets with isolated leptons, 
plus the total event rate, should allow for an extraction of the gluino mass. 

The soft component of signal will come from XiXi ^ X1X2 ™d X^Xi production. These 
events will be followed by 3-body decays to various final states, but since the visible com- 
ponents of the signal are much softer than that from gluino pair production, these events 
will be harder to see above SM background levels. With judicious cuts, the soft component 



might also be visible at some level (e.g. xfx2 ^ 3^ + i?™*^*^) [46]. 
Point B 

Point B is the same as point A, except that in this case the gaugino mass Mi has been 
raised to 195 GeV so that the xf — x? mass gap shrinks to only 13 GeV. Since n is quite 
large, the Xi remains nearly pure bino-like, but the relic density problem is solved via 
bino-wino co-annihilation. This case will again give a hard component to the LHC new 
physics signal from gluino pair production, but this time the m{i^i~) distribution will 
have an edge only at 13 GeV. When compared to any gluino mass reconstructions, this 
would indicate a violation of gaugino mass unification at the GUT scale. In addition, the 
small X2 ~ Xi mass gap suppresses 3-body decays such as X2 ~^ XiQQ and Xi^^ relative to 
any kinematically-allowed 2-body decays such as the loop-induced process X2 ~^ Xil 01 • 
Thus, the radiative X2 decay to photon X2 ~^ Xi7 can become large |24|: in this case, it 
reaches 10%. The final state 7 will be somewhat soft if the X2 velocity. But if X2 

is moving fast as a result of production from cascade decays, then hard, isolated photons 
should occasionally be present in the SUSY collider events. 

Point C 

The point C parameters are listed in Table In this case, mi6(l,2) has been lowered 
far below mi6(3) so that the first two generations of scalars are degenerate, but with a 
lower mass than third generation scalars. The Higgs mass splitting leads to a large RGE 
5-term, which drives ur and cr to very low masses ~ 98.3 GeV. This is at the edge of 
LEP2 exclusion. The relic density problem is solved because XiXi ~^ QQ via qr exchange 
and neutralino-squark co-annihilation act to reduce the relic density. The cross section for 
production of two flavors of extremely light squarks is extremely large at LHC. Normally 
one would expect characteristic dijet-|-£^™'**** events since qr — > qxi- However, in this case 
the mass gap m^^ — m^o ~ 18 GeV, so both the jets and E!^^^^ will be very soft. Gluinos 
and other squarks will also be produced at large rates, although the g uur, ccr decays 
are dominant. While left-squarks may decay with large rates to X2 Xi^) we note that 
X2 — ^ UUR and ccr is also large, leading again to relatively soft jet activity. In spite of the 
soft jet activity, the scenario should be easily seen at LHC, since ql q'Xi occurs at a 
large rate, and Xi ~^ cfeX? occurs at 43% branching fraction (enhanced by the relatively 
light left-sleptons). This can lead to a large same-sign (SS) dilepton rate from pp ^ ulul 
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parameter 


A 


B 


C 


D 


E 




9202.9 


9202.9 


5018.8 


2976.5 


5877.3 


"^1/2 


62.5 


62.5 


160 


107.0 


113.6 


^0 


-19964.5 


-19964.5 


-10624.2 


-6060.3 


-12052.6 


mio 


10966.1 


10966.1 


6082.1 


3787.9 





tan/3 


49.1 


49.1 


47.8 


49.05 


47.4 


Md 


3504.4 


3504.4 


1530.1 


1020.8 


— 


Ml 


— 


195 


— 


— 


— 


mi6(l,2) 


— 


— 


603.8 


— 


— 


/t 


0.51 


0.51 


0.49 


0.48 


0.49 


fb 


0.51 


0.51 


0.41 


0.47 


0.49 


fr 


0.52 


0.52 


0.47 


0.52 


0.49 


M 


4179.8 


4186.3 


1882.6 


331.0 


865.3 


mg 


395.6 


395.4 


495.5 


387.7 


466.6 




9185.4 


9185.4 


622.1 


2970.8 


5863.0 




9104.1 


9104.2 


98.3 


2951.4 


5819.2 




zolo.i 


zdlU.o 


11)48.4 


AO A C 

434.5 


C\A A 'y 

944.7 




z / zo.i 


z 1 14. y 


ioy4.u 


o4y.o 


1 A c;o 1 
i4oz. / 




yioi.y 


yioz.U 


oil.y 


OOP; ^ Q 


OOOO.O 




c/OZO. 1 






OUUa.U 


Oa'-±'J.O 


m-± 


128.8 


128.8 


165.7 


105.7 


141.3 




128.6 


128.1 


165.1 


105.1 


140.9 




55.6 


115.9 


80.2 


52.6 


65.7 


rriA 


3273.6 


3266.0 


1939.9 


776.8 


177.8 




125.4 


125.4 


123.2 


111.1 


113.4 


Al 

BF{b sj) 

BF{Bs ^ 

(^sciXiP) [pb] 


423 

220 

3.0x10-" 

3.3x10-4 
5.0x10-12 

5.1x10-12 

5.0x10-9 

4.4x10-9 

1.3 X 10-15 


0.09 
0.08 

3.0x10-4 

3.3x10-4 
5.0x10-12 
5.0x10-12 
5.0x10-9 
4.4x10-9 

1.9 X 10-1^ 


0.11 
0.11 

6.2x10-4 

3.7x10-4 
3.0x10-1" 
2.8x10-1" 
11.8x10-9 

6.9x10-9 

1.5 X 10-6 


0.10 
0.06 

1.9x10-4 

4.0x10-4 
2.2x10-1" 
2.2x10-1" 
5.8x10-8 
6.2x10-8 

2.7 X 10-9 


0.15 
0.08 

2.5x10-4 

2.2x10-4 
4.1x10-11 
4.1x10-11 
2.0x10-5 
2.0x10-5 

5.3 X 10-8 



Table 1: Masses and parameters in GeV units for five benchmark Yukawa unified points using 
Isajet 7.75 and mt = 171.0 GeV. The upper entry for the ^^l^^ come from IsaReD/Isatools, 
while the lower entry comes from micrOMEGAs; (t(x5p) is computed with Isatools. 



production, along with a large asymmetry in ++ SS dileptons over SS dileptons (which 

occur from d^dL production). This scenario may also be subject to exclusion by analysis 
of Fermilab Tevatron data. We further note that point C is naively excluded by direct 
dark matter search limits. These latter limits depend on an assumed standard local relic 
density mass and velocity distribution, so that the limits can be avoided if one postulates 
that we live in a local underdensity of dark matter. 
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Point D 



Point D is an example of a compromise solution, where we allow mig as low as 3 TeV at 
some expense to Yukawa unification (here, Yukawa unification is good to only ~ 10%) in 
order to allow for neutralino annihilation through the light Higgs /i-resonance (neutralinos 
can still annihilate through the light h resonance for higher uiiq values; it is just that the 
relic density can't be pushed as low as ^j^ah'^ ~ 0.1). This scenario is extremely predictive, 
with gluinos around 350-450 GeV, so again we expect LHC events to be dominated by 
gluino pair production. As in the case of point A, the gg events will be followed by 3-body 
decays to 6-jet rich final states. A dilepton mass edge at m^o — m^o ~ 53 GeV should be 
visible since X2 ~^ Xi^^e" at 3.3% branching fraction. The ii weighs only 434 GeV in this 
case, and bi is at 849 GeV, so it may be possible to detect some third generation squark 
pair production events. The top squark decays to bxf with a 50% branching fraction, and 
also has significant branching fractions to txi, a-nd &X2^ final states. The 61 dominantly 
decays to bg and Wti final states. Moreover, the heavy Higgs bosons A^, H° and have 



masses around 780 GeV and should be detectable at LHC ||48|| . 
Point E 

Point E is a Yukawa-unified solution that solves the DM abundance problem via neutralino 
annihilation through a 178 GeV pseudoscalar A resonance. The combination of light A 
and large tan /3 leads to a branching fraction Bg fJ-^fJ-^ which is excluded by recent 
CDF analyses. If we are allowed to somehow ignore this (possibly via other fiavor-violating 
interactions), then the scenario would be at the edge of observability via Tevatron searches 
for A, H ^ T^T~ and bb, which at present exclude tua ^ 170 GeV [|^]. The LHC (and 
possibly soon also the Tevatron) would easily see the rather light spectrum of Higgs bosons. 
Gluinos can be somewhat heavier in this case compared to points A and D, ranging to over 
a TeV. However, in point E as listed, with a 467 GeV gluino, the gluino pair production 
signatures will be rather similar to those of point A: rich in 6-jets, with a visible dilepton 
mass edge at 75 GeV. 



5. Summary and conclusions 

In this paper, we have presented a number of new results. 

1. First, we verified former results presented in Ref. |l^] that Yukawa unified models can 
be generated with updated Isajet spectra code and an updated value of the top quark 
mass nit = 171 GeV. Using both random scans and the more efficient MCMC scans, 
we find that models with excellent Yukawa coupling unification can be generated in 
the HS model if scalar masses are in the multi-TeV range, while gaugino masses are 
quite light, and the Xi is slightly above the current LEP2 limit. The models require 
the Bagger et al. boundary conditions if ^ > such that A^ = ^rn^Q = 4m|g, and 
^0 < in our convention. The spectra generated is characterized by three mass scales: 
multi-TeV first and second generation matter scalars, TeV scale third generation and 
Higgs scalars and 100-200 GeV light charginos and gluinos of order 350-450 GeV. 
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The relic density is typically 30-30,000 times above the WMAP measured value. As 
a solution, we propose i). hypothesizing an unstable neutralino x5 which decays to 
axino plus photon, ii). raising the GUT scale gaugino mass Mi so that bine- wine 
co-annihilation reduces the relic density or in), lowering the first/second generation 
scalar masses relative to the third so that neutralinos can annihilate via light qa 
exchange and neutralino-squark co-annihilation. We regard the first of these solutions 
as the most attractive, and the third is actually susceptable to possible exclusion by 
analyses of Fermilab Tevatron signals in the case of just two light squarks. 

2. Using an MCMC analysis, we find a new class of solutions with mie ~ 3 TeV, where 
neutralinos annihilate through the light higgs h resonance. This low a value of miQ 
typically leads to Yukawa unification at the 5-10% level at best. 

3. We find we are able to generate solutions with low and low iua as did the BDR 

group. The solutions generated by the Isajct code with low /x, low niA and uiiq ~ 3 
TcV tend to have Yukawa unification in the 20% range or greater. We were able to 
generate a class of solutions with excellent Yukawa unification and rriiQ ranging up to 
6 TeV, where the DM problem is solved by neutralino annihilation through a 150-250 
GeV A resonance. The combination of large tan /? and low ruA gives a Bg ^ 
branching fraction at levels beyond those allowed by the CDF collaboration. 

We also present a Table of five benchmark solutions suitable for event generation, and 
for examination of collider signals expected at the LHC from DM-allowed Yukawa-unified 
SUSY models. Based on this work, we are able to make several predictions, if the Yukawa- 
unified MSSM is the correct effective field theory between Mqut and Myjgak- We would 
expect the following: 

• New physics events at the CERN LHC to be dominated by gluino pair production 
with nig ~ 350-450 GeV. Since tan/3 is large, the final states are rich in 6-jets, 
and the OS/SF isolated dilepton invariant mass distribution should have a visible 
edge at m^o — m^o ~ 50-75 GeV because the X2 always decays via 3-body modes. 
Squarks and sleptons are likely to be very heavy, and may decouple from LHC physics 
signatures. 

• We would predict in this scenario that the {g — 2)^ anomaly is false, since in Yukawa- 
unified SUSY models with large uiiq, the SUSY contribution to the muon QED vertex 
is always highly suppressed. 

• While SUSY should be easily visible at the LHC for Yukawa unified models, we 
would predict a dearth of direct and indirect dark matter detection signals. This 

is because the typically large values of n and scalar masses tend to suppress such 
signals. However, in the CDF-excluded case of point E, the direct and indirect DM 
signals may be observable. Points C and D also have low but observable direct DM 
detection rates, since scalars are not too heavy. 
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